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Abstract 


Application  of  superconducting  electrical  machinery  to 
ship  propulsion  yields  light  weight  and  small  volume  propul- 
sion plants  with  increased  system  efficiency,  reliability 
and  lower  cost.   Such  propulsion  systems  might  be  the  only 
feasible  answer  for  high  performance  ship  types. 

A  design  method  for  two-pole  synchronous  machine  with 
superconducting  field  winding  is  developed  which  will  survive 
even  under  the  worse  fault  conditions.   Iron  shielded  machines 
with  the  effects  of  higher  harmonics  and  partial  winding 
angles  are  considered.   The  method  is  then  adapted  to  a 
computer  program  for  ease  of  obtaining  minimum  volume  and 
weight  machines.   The  results  are  compared  with  20,000  HP 
conventional  propulsion  generator.   It  is  shown  that  great 
weight  and  volume  savings  are  possible. 
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NQMENCLATURE 

Subscripts 

a  armature  winding 

f  field  winding 

r  radial 

6  azlrauthal 

i  inside  (for  field  quantities — at  a  radius  less 
than  inner  radius) 

i  inner  (when  used  with  a  dimension) 

o  outside  (when  used  with  a  field  quantity) 

o  outer  (when  used  with  a  dimension) 

s  field  quantities,  generated  by  a  sheet  current 

a  phase  a 

b  phase  b 

c  phase  c 

d  direct  axis-armature 

q  quadrature  axis-armature 

sd  direct  axis-damper  shield 

sq  quadrature  axis-damper  shield 

Symbols 

B'  magnetic  flux  density 

^sat  saturation  flux  density  of  shield  material 


Bj^^^2     mean-squared  field  experienced  by  armature 
conductors 


av 


w       diameter  of  primary  conductor 
dg^      dissipation  density  of  shield  material  (watts/lb) 


E-:,  open  circuit  voltage  (RMS) 

E„"  voltage  behind  subtransient  reactance 

e^"  voltage  behind  subtransient  reactance  (per  unit) 

H   '  magnetic  field  density 

i  current 

I  rated  armature  current  (RMS) 
a 

J  current  density 

J  rated  armature-current  density 
a 

J  rated  field-current  density 

h  shield  structure  thickness 

h,  primary  shield  thickness 

i  straight-section  length 

i  total  machine  length 

I  effective  length  for  calculating  armature 
self  inductance 

Ir,  effective  length  for  calculating  field-winding 
self -inductance 

i  effective  length  for  calculating  field-to- 
m  ^^  >— ' 

armature  mutual  inductance 


Z^  effective  length  for  calculating  damper  self- 

inductance 


'd 


a    ,      effective  length  for  calculating  damper-to- 
armature  mutual  inductance 

A-,      effective  length  for  calculating  field-to- 
damper  mutual  inductance 

a  length  of  armature  conductors 

i    .  effective  length  for  calculating  eddy-current 

losses 

i^Yi  effective  length  for  calculating  shield  losses 

L        self  inductance,  armature  phase 
^  8 


L«  self  inductance,  field 

L  ^  mutual  inductance,  phase-to-phase 
ab 

L  ^  mutual  inductance,  armature-to-damper 
ad 

L»,  mutual  inductance,  field-to-damper 

M  magnitude  of  mutual  inductance,  field-to-armature 

M  mutual  inductance,  field-to-phase-a 
a 

M  magnitude  of  mutual  inductance,  damper-to-pha£'.e-a 

N  ,  number  of  turns  in  armature  winding 
at 

N-^,  number  of  turns  in  field  winding 

p  number  of  pole  pairs 

P  rated  volt-amperes 

P  conduction  losses 
c 

P  ,  eddy-current  losses 
ed 

P  ,  shield  dissipation 
sh  ^ 

R  inner  radius,  shield 
s  ' 

R  outer  radius,  shield 
OS  ' 

R,  outer  radius,  damper  shield 

V,  rated  terminal  voltage 

X„  synchronous  impedance  (in  ohms) 
a 

X  per-unit  synchronous  impedance,  normalized  to 
internal  voltage,  E^ 

X,  per-unit  synchronous  impedance,  with  V,  as 
base  voltage 

X,'  per-unit  transient  reactance  (direct  axis) 

X , "  per-unit  subrtansient  reactance 
d 

y   =   Rf^/Rj,^ 


5.  skin  depth 

0  angular  displacement 

0  -  field-winding  included  angle 
wf  .  ^  ^ 

0  ^   =   pS  ^ 
wfe       wf 

0  armature  winding  angle 
wa  &   & 

0     =  p0 
wae       wa 

a  magnitude  of  maximum  radial  stress  under  three 
phase  fault  from  load 

a  radial  stresses 

T  sheer  stresses 

p  density 

a  conductivity  of  armature  conductors 

a^  conductivity  of  damper-shield  material 

<p  rotor-axis  displacement  from  axis  of  phase-a 

0)  angular  frequency,  also  used  for  deflections 

u)  electrical  angular  frequency  =  pw 
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INTRODUCTION 

Advantages  of  electrical  ship  propulsion  systems  were 
realized  years  ago.   Some  of  them  are:   arrangement  flexib- 
ility, ease  of  automation  and  possible  integration  to  ship 
service  power.   Major  drawbacks  were  the  relatively  high 
cost  and  weight  of  system  components. 

Over  the  past  ten  years  or  so  many  studies  have  been 
done  on  the  possible  applications  of  superconductors  to 
electrical  power  equipment,  with  most  effort  devoted  to 
transformers  with  decreasing  attention  given  to  transmission, 
ac  machines,  and  dc  machines  respectively.   The  development 
of  high-field  superconducting  wire,  suitable  for  fabrication 
into  windings  for  electric  machines,  high  vacuum  seals  and 
rotating  dewar  technology,  made  applications  of  supercon- 
ductors possible  for  rotating  electric  machines. 

Among  the  many  different  types  of  superconductors , the 
most  useful  one  is  type  II,  because  of  its  ability  to  carry 
a  high  current  density  in  a  high  magnetic  field,  and  the 
first  major  consequence  of  this  upon  machine  design  is  that 
we  may  break  free  from  the  magnetic  limitations  of  iron. 
Elimination  of  iron  will  reduce  the  machine  weight  and 
provide  more  armature  conductors  for  a  given  space,  thus 
increasing  the  power  rating  per  unit  volume.   If  the  flux 
density  at  the  armature  conductors  is  increased  by  a  factor 
of  3  over  that  in  a  conventional  machine,  the  generated  emf 
per  unit  length  of  conductor  is  thus  increased  3  times,  if 
the  armature  diameter  and  the  rotor  speed  are  left  unchanged. 
Therefore,  for  the  same  output  power  and  current,  the 
machine  length  can  be  reduced  by  the  factor  3s   with  a 
corresponding  reduction  in  machine  volume  and  weight. 

Application  of  superconductors  to  rotating  machines  is 
not  without  problems.   First  of  all,  they  will  operate  only 
at  extremely  low  temperatures.   Secondly,  perhaps  the  most 
important  one  is  that  time  varying  magnetic  fields  produce 

11 


losses  v/hich  are  prohibitive  at  power  frequencies.   These 
losses  are  the  limiting  factor  for  the  application  of  super- 
conductors to  field  windings.   Therefore,  the  configuration 
of  a  superconducting  machine  must  be  a  superconducting 
stationary  or  rotating  field  winding  and  a  normally  conduct- 
ing ac  armature  winding,  which  may  also  be  rotating  or 
stationary.   Configuration  for  rotating  armature  has  inherent 
problem  of  transmitting  power  across  sliding  contacts. 
With  the  advent  of  rotatable  liquid  helium  transfer  coupling 
a  rotating  superconducting  field  winding  became  practical, 
thus  eliminating  the  necesity  for  rotating  the  armature. 

Different  designs  and  configurations  have  been  studied  in 

567 
different  countries,  among  them  synchronous  machines  '  '  , 


1 
3 


•5  -^   l^ 

dc  machines  ,  homopolar  machines  '  ,  double  armature  ma- 

1  2 

chines  .  and  the  toroidal  reciprocating  alternator  .   Some 

of  them  are  already  built  or  under  development. 

Application  of  superconducting  electric  machinery  to 
ship  propulsion  systems  has  been  studied  over  approximately 
the  past  seven  years  '  '  '   '   .   A  superconducting  propul- 
sion system  is  shown  in  Figure  1.   The  required  set  of 
performance  characteristics  from  a  propulsion  plant  will 
depend  mainly  on  the  ship  type  being  considered.   Such 
characteristics  are  tied  with  the  options,  control  and  type 
of  electric  machinery.   For  example,  bulk  carriers  which 
run  at  their  full  speed  Q5%   of  their  operating  time  will 
not  require  speed  control  as  much  as  a  naval  combatant 
vessel,  and  the  option  box  may  be  omitted  and  a  simple 
synchronous-synchronous  drive  in  which  the  propeller  speed 
can  be  controlled  by  varying  prime  mover  speed,  may  be 
used.   Destroyers,  which  operate  over  their  full  speed  range, 
will  require  more  sophisticated  controls  and  the  option  box 
may  be  filled  with  a  cyloconverter  if  a  synchronous  system 
is  used  or  dc  machinery  might  be  more  appropriate.   It 
should  be  pointed  out  that  the  choice  of  thruster  and  prime 
mover  will  also  affect  the  electric  machinery. 

12 
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Recently  some  studies  have  been  done  on  conventional 

12  13  l4 
electric  propulsion  of  bulk  carriers   »  -^ »   .   Sixteen 

tankers,  container  ships  and  LNG  carriers  are  being  built 
or  are  under  design.   Conventional  electric  propulsion 
systems  have  been  selected  due  to  the  fact  that  flexibility 
of  propulsion  systems  may  lead  to  possible  savings  in  the 
volume  of  the  plant.   Gas  turbine  generator  sets  can  be 
placed  on  the  main  deck,  reducing  the  inlet  and  exhaust 
ducting  volume  and  flow  pressur-e  drops,  thus  improving  gas 
turbine  performance.   Added  to  these  advantages,  if  a  super- 
conducting drive  is  used  even  more  savings  are  possible. 

15 
It  can  be  shown  that    the  reduction  of  specific  machinery 

weight  from  30  to  10  Ibs/SHP  will  permit  an  increase  in 

maximum  speed  from  35  to  50  knots  with  no  sacrifice  in 

payload  or  endurance  for  a  destroyer  type  ship.   Introduction 

of  gas  turbines  has  made  specific  weights  somewhat  lower, 

but  introduced  new  problem  areas.   Newly  developed  ship 

types  such  as  hydrofoils,  SES  and  catamarans  which  require 

unusual  propulsion  plant  arrangements,  suffer  greatly  with 

conventional  propulsion  plants.   With  increased  power  levels 

required  by  high  performance  ship  types,  superconducting 

electric  propulsion  plants  may  be  the  only  feasible  answer. 

The  purpose  of  this  thesis  is  to  obtain  a  design  method  for 

synchronous  machines  which  can  be  used  to  realize  possible 

advantages. 
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Chapter  1 
DEVELOPMENT  OF  A  MODEL 

A  Model  Configuration 

The  geometrical  model  used  for  this  analysis  is 

illustrated  in  Figure  2  and  Figure  3,    with  approximate 

account  taken  of  the  effects  of  end  turns  as  indicated  by 

5 
measurements  made  on  the  first  MIT  experimental  machine  . 

The  axial  length  of  one  set  of  armature  end  turns  is  taken 

to  be  equal  to  the  average  armature  winding  diameter  in  the 

active  section  divided  by  the  number  of  pole  pairs,  and  the 

radial  thickness  of  the  end  turns  is  assumed  equal  to  twice 

the  radial  thickness  of  the  armature  winding  in  the  active 

region.   These  relatively  large  end  turns  are  assumed 

because  there  is  more  copper  in  the  active  region  and  more 

volume  is  needed  in  which  to  cross  the  phase  conductors 

over  each  other. 

The  machine  has  a  smooth  outer  shell  made  of  either 
laminated  iron  or  a  highly  conducting  material,  such  as 
copper  or  aluminum,  in  order  to  provide  a  uniform  boundary 
condition  and  confine  the  magnetic  field  within  the  machine 
In  the  case  of  conductive  image  shield,  circulating  image 
currents  serve  to  confine  the  magnetic  field. 

There  is  a  thin  highly  conductive  shield  between  the 
field  winding  and  armature  which  rotates  with  the  field 
winding.   The  purpose  of  this  shell  is  to  intercept 
asynchronous  magnetic  fields  due  to  load  imbalance  and 
space  harmonics  of  the  armature  field  and  take  large  fault 

-1   -1  o 

torques  and  radial  loads  under  fault  conditions  '   .   It 
may  also  intercept  thermal  radiation  from  the  room  tempera- 
ture armature  or  a  separate  thermal  radiation  shield  may  be 
used.   In  any  case,  the  essential  purpose  of  this  shell  is 
to  protect  the  field  windings.   Since  this  shell  will  not 
be  strong  enough  to  take  large  fault  loadings,  it  should 
possibly  be  supported  by  a  shell  made  out  of  a  stronger 
material.   In  this  work,  the  load  carrying  structure  will 

15 
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be  called  the  shield  structure  and  the  shell  v:hich  provides 
actual  shielding  will  be  called  the  primary  shield.   Shown 
in  Figures  2  and  3,  the  primary  shield  is  inside  the  shield 
structure,  but  it  may  very  well  be  outside  the  shield 
structure.   Consequences  of  these  configurations  will  be 
discussed  in  a  later  section. 

The  electrical  configuration  is  shown  in  Figure  4. 
There,  are  a  total  of  six  windings,  three  phase,  one  field 
and  a  direct  and  a  quadrature  axis  damper  winding.   The 
electrical  parameters  of  the  machine  are  the  coefficients 
of  the  flux-current  relationship  which  can  be  expressed  in 
matrix  form   '   : 
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Figure  4 

Electrical  Configuration 
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This  may  be  expressed  by: 


A 
a 

^b 

X 
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= 

K& 

/sq. 
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'ab 


'ab 


'ab 
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M   cosp(|)        M   cosp((t)-27T/3)         M   cosp  (({)+27t/3  ) 


M      cosp(J)      M      cosp((|)-2iT/3)      M      cosp  ((J)+2tt/3  ) 

So  o. 


M      sinpcj)      M      sinp((})-27T/3)      M      sinp  ((j)+2Tr/3) 
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By  doing  so,  it  is  assumed  that  all  rotor-to-armature 
mutual  inductances  vary  sinusoidally  with  rotor  position. 
This  is  equivalent  to  neglecting  higher  order  harmonics. 
It  should  be  emphasized  that  higher  order  harmonics  in 
phase-to-phase  mutual  inductances  were  not  neglected  in  this 
study.  20 


IB  Field  Analysis 

?or  the  purposes  of  this  study,  a  two  dimensional 
representation  is  used.   The  effects  of  finite  length  and 
end  turn  geometry  will  be  included  in  an  approximate  fashion 
In  the  straight,  active  section  of  the  machine,  currents 
are  parallel  to  the  axis  of  the  machine  and  effective 
armature  and  field  current  densities  are  related  to  terminal 
current  by 


2  ^a  ^at 

^   e  ^  R  '^d-x'^) 

wae  ao 


2  ^f  ^ft 
Jf.  =  ^  p^^   p  (2) 

Ve  \o^l-y^ 


The  field  and  armature  winding  are  assumed  to  have  a 
uniform  density  between  their  inner  and  outer  radii,  and 
to  fill  the  wedged  shape  space  within  a  winding  angle. 
Expressions  for  the  magnetic  fields  generated  by  these 
currents  are  derived  in  references  l6  and  17,  and  are 
summarized  in  Table  1.   Wherever  double  signs  (±  or  +) 
appear,  the  upper  sign  refers  to  the  laminated  iron  shield 
case,  the  lower  to  the  conductive  image  shield  case.   Unless 
otherwise  specified,  all  current  densities  take  space  factors 
into  account. 
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TABLE  1 


MAGNETIC  FIELDS 
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Table  1 

Magnetic  Fields  (continued) 
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Tab]e  1 

Magnetic  Fields  (continued) 
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Inductance  Parameters  • 

Self  and  mutual  inductances  of  and  between  the  various 
windings  are  calculated  by  integrating  flux  linkage  of  the 
fields  given  in  Table  1,  over  the  area  of  a  winding. 
Reference  l6  gives  such  a  calculation  and  the  results  taken 
from  this  reference  are  summarized  in  Table  2. 

The  electrothermal  shield  has  been  treated  as  a 
fourier  series  of  sinusoidally  wound  damper  windings.   The 
justification  for  this  is  that  each  space  harmonic  field 
will  excite  currents  in  this  shell  that  are  sinusoidally 
distributed  in  angle.   Further  more,  each  space  harmonic  is 
orthogonal  to  all  but  one  of  the  ficticious  windings.   Each 
winding  is  characterized  by  a  self  inductance  and  mutual 
inductances  with  field  and  armature  windings.   The  mutual 
inductance  between  any  two  of  these  damper  windings  is  zero. 
Included  in  the  expressions,  Nnp  is  fictitious  number  of 
turns  of  damper  winding. 

The  expressions  given  in  Table  2  have  been  left  with  the 
effective  length  noted  differently  for  each  parameter. 
These  lengths  may  not  correspond  to  any  physical  lengths  of 
machine  parts,  because  of  the  effects  of  end  turns.   Further 
more,  since  end  turns  of  different  coils  in  a  machine  will 
have  different  lengths,  we  miay  expect  the  coils  to  have 
differing  effective  lengths.   Certain  assumptions  utilized 
when  evaluating  effective  lengths  are  discussed  in  section 
D  of  this  chapter. 
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. TABLE  2 
INDUCTANCES 
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ID   Machine  Rating  Equations 

The  volt-ampere  rating  of  a  machine  is  given  by 


P  =  3V,I^  =  3E^I^(g^) 


(3) 


where  V ,  ,  E„,I  ,  are  RMS  quantities 
The  internal  voltage  E^  is  given  by 
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U^^MI^ 
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(4) 


V. 
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E^'   Kl-x^^cos^y^   -  X  sin^ 
r        a  a 


(see  Fig.  5)  (5) 
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a 


Figure  5 

Voltage-current  relationship 
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where  x  is  synchronous  reactance  normalized  to  internal 
voltage 
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They  may  be  combined  to  give 


P  =  2i  ^    1   J   J  (i_yP+2)R   P+2    2-p 
75-   e  0  m  f  a^   -^    ^  fo    ao 


fe      ^ 

wae 
2 

sin 

^wfe 
2 

'ip 


where  C,   is  given  from  Table  2 
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sin 


Ip 


where  C  is  given  by 

ne 


sm 


C  =    E 
n  odd 


wae 


(1-cos 


2n7r 


LU 


n 


snp 


and  C    can  be  obtained  from  Table  2. 
snp 

Transient  reactances  are  calculated  with  the  assumption 
that; 

a)  Field-winding  flux  is  constant 

b)  Damper-winding  current  is  zero 


These  reactances  are 


16,17 


'^'  ^  ^^^'"2'(v4b^^f^ 


(8) 


q     d 


where  X,  is  per  unit  synchronous  reactance 


■y  _     B.       a   _   y  1 


(9) 


Subtransient  reactances  are  calculated  assuming  that 

flux  linked  by  the  damper  shield  is  constant.   These  are 

u   16,17 
given  by   ' 


X  ••  =  X  ''  =  X,(l-^ 


3 


M  2 

^vC^^di 


(10) 


Since  the  shell  is  symmetric,  direct  and  quadrature  axis 
subtransient  reactances  are  the  same. 
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IE  Machine  Losses 

There  are  mainly  three  different  losses  in  synchronous 

16 
machines   : 


a)   Current-conduction  losses: 
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b)   Eddy-current  losses: 
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where  B    is  the  mean-squared  magnetic  field  seen  by  the 
av 

armature  conductors: 
2      1 
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ao  ai 


J     B^.  ''.dr 
ai 


(13) 


where  B^„  is  given  in  Table  1.   If  the  first  harmonic  of 
B^  is  considered,  for  a  two-pole  (one  pole  pair)  iron 
shielded  machine: 
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o)   Shield  losses: 

If  an  image  shield  is  used,  circulating  currents 
will  produce  losses: 

P  ,  =  — ®^ ^— ^  (15) 

sh 

OSS 


where  B^  is  B.  evaluated  at  shield  inside  radius,  and  6 

QS        ©  '         S 

is  the  skin  depth  of  the  shield  material. 

If  an  iron  shield  is  used,  eddy-currents  and  hysterisis 
will  produce  losses: 

P  ,   =  7r(R   ^-R  ^)9  V.  P  V,  d  ,  '  ..rs 

sh      OS    s  ^sh   sh   sh  (16) 


where  d  ,  is  shield  dissipation  per  unit  mass  of  shield 
sh  j^      t- 

material.   The  effective  lengths  will  be  discussed  in 
Appendix  B. 
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IF  Fault  Effects 

Although  elimination  of  iron  reduces  the  machine  vol- 
ume and  weight,  it  degrades  the  transient  behaviour  of 
superconducting  machines.   Large  fault  torques  and  stresses 
will  be  imposed  upon  the  field  winding.   In  addition  to 
fault  loading,  there  will  be  heat  dissipation  and  field 
current  rise,  which  may  drive  superconducting  field  winding 
normal.   The  protection  needed  will  be  provided  by  a  suitable 
shield  around  the  field  winding.   The  expressions  for  fault 
loadings  under  various  fault  conditions  are  derived  in 
references  1  and  18.   For  us  the  maximum  values  are  import- 
ant in  order  to  design  a  machine  which  can  withstand  such 
loads. 

The  maximum  radial  loads  imposed  upon  the  primary 

shield  will  be  under  three  phase  short  circuit  from  load. 

1  1  ft 
These  are  given  by  '   : 
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(17) 


where  H  ^   and  H  .  are  the  fields  outside  and  inside  the 
shell  respectively  after  the  fault: 
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and  H   is  the  armature  field  at  primary  shield  radius: 
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if  P*2 


and  H^  is  the  field  generated  by  field  windings  at  primary 
shield  radius: 
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It  should  be  noted  that  only  the  first  harmonic  fields 
have  been  taken  into  account. 

When  equations  (18),  (19) j  and  (20)  are  combined  with 
(17)  and  the  maximum  stress  evaluated  when  wt„='n'. 


aj,=  2  y^H^d+f )  [H^  (1+f  )+H^  qIcos^  (p  6-  cf)^) 


(21) 


The  maximum  torque  will  be  applied  on  the  primary 

shield  under  line-to-line  short  circuit  from  load.   This 

1  1  8 
torque  is  given  by  '   : 


■  e  "2 
f 
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(22) 


and   the   maximum  value   will    occur   at    wtp.=  ■n/2,for   which 


,,2 
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In  above  derivations  a  factor  e„",  voltage  behind  the 
subtransient  reactance,  was  introduced  to  correct  fault 
loads  from  open  circuit  to  load  in  an  approximate  fashion. 
Figure  6  shows  a  method  to  calculate  e„"  from  which 


F  " 


=  Xl+X^"  sinl//)^  +  (X^"  .  cosy/)^     (24) 


j(Xa-Xa)Ia 
J  ^a  la 


Figure  6 
Vector  Diagram  of  Internal,  Terminal  and  Reaction 

Voltages 
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Since,  primary  shield  will  not  be  strong  enough  to 
resist  fault  loadings,  it  must  be  supported  by  a  shield 
structure.   Stresses  induced  under  radial  and  torque  loads 
are  derived  in  Appendix  A, 

Under  fault  conditions,  the  superconducting  field 
winding  may  be  driven  normal,  due  to  the  heat  dissipation 

and  the  field  current  rise.   Perhaps  the  most  important  one 

1  R 
is  the  field  current  rise    in  both  dc  and  ac  levels.   Allow- 
able field  current  rise  will  depend  on  the  rated  operating 
point  with  respect  to  superconductor  short  sample  curve. 
Figure  7  illustrates  the  allowable  field  current  rise  for 

a  given  operating  point.   Ac  component  of  the  field  current 

1  R 
rise  will  be  small    compared  to  the  dc  component  and  may 

be  neglected.   Field  current  will  rise  exponentially  and 

its  final  value  will  depend  on  primary  shield  time  constant, 

1  1  R 
machine  reactances  and  circuit  breaker  opening  time  '   : 


1,  =  i,„[X  .  ^  (,.e-t/T,",, 


(25) 


where  i^  and  i„   are  the  instant  and  initial  values  of 
1       fo 

field  current  respectively;  and,  T . "  is  related  to  the 


primary  shield  time  constant  by 
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and  for  a  thin  primary  shield,  shield  time  constant 
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^  -y  a  t   R,  [1± 
SI    2p^o    s   t 


fRj2p 
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] 


(26)  . 
18,16 

(27) 


This  time  constant  will  also  determine  the  machine  swings 
under  transient  conditions,  but  this  was  not  included  in 
this  study. 
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Figure  7 
Allowable  Field  Current  Rise 
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IG   Design  Constraints 

A  superconducting  coil  is  characterized  by  its  short 
sample  test  curve.   At  a  given  operating  temperature,  the 
relation  between  flux  density  and  current  density  which  can 
be  carried  without  driving  superconductor  normal,  constitutes 
such  a  curve.   Any  chosen  operating  point  should  be  to  the 
left  of  this  curve.   In  order  to  provide  stability  against 
disturbances  of  electrical,  thermal,  or  mechanical  nature 
which  may  occur,  a  material  of  high  electrical  and  thermal 
conductivity  is  included  in  the  conductor.   Copper  is 
usually  used  for  this  purpose.   High-purity  aluminum  has  not 
been  accepted  because  of  the  problem  of  making  a  metal- 
lurgical bond  with  the  superconductor  although  it  has  higher 

electrical  conductivity.   Recently  a  metallurgical  bond  has 

22 
been  made    but  the  brittle  nature  of  it  still  limits 

aluminum  as  a  stabilizer.   Figure  8  shows  Nb  Ti  superconductor 
short  sample  test  data,  with  the  state  of  the  art  line  given 
in  reference  21.   To  the  right  of  this  line  small  disturb- 
ance may  drive  superconductor  to  normal,  according  to  the 
stability  criteria  given  in  this  reference.   Although  such 
a  line  is  quite  arbitrary,  it  defines  a  region  which  is 
experienced  with  other  magnets.   It  should  be  pointed  out 
that  the  higher  flux  densities  are  possible  with  Nb-.  Sn 
conductors,  but  they  are  more  expensive  than  Nb  Ti 
conductors  and  cannot  be  fabricated  into  complex  windings. 

Another  major  limitation  comes  from  the  saturation 
considerations  of  the  iron  shield  used  to  confine  magnetic 
fields  inside  the  machine  volume.   The  shield  may  be  locally 
saturated  or  the  maximum  flux  density  carried  by  the  shield 
may  exceed  saturation  flux  density  of  the  iron.   These 
considerations  will  be  discussed  in  detail  in  Appendix  B 
when  the  computer  program  is  described. 

Armature  current  density  will  be  determined  by  the 
losses  and  cooling  method  employed.   A  figure  of  3«0  x' 10 
amp/m  could  be  accepted  for  a  nitrogen  cooled  armature. 
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Figure  8 
Superconductor  Short  Sample  Test  Data 
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structural  integrity  of  the  machine  structure  will  also 
constitute  a  problem  especially  when  it  is  degraded  under 
fault  conditions.   The  most  critical  part  of  the  machine 
will  be  the  shield  structure  and  must  be  designed  to  take 
large  fault  loadings.   This  will  be  discussed  in  Appendix 
A  and  B,  in  detail. 
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Chapter  2 
RESULTS 

2 A   Parameters  Used  in  the  Design  Study 

For  the  purpose  of  this  study,  calculations  were 
performed  for  a  20,000  HP,  two-pole  synchronous  generator. 
Input  parameters  to  the  computer  program  are  given  in  Table 
3.   First  harmonic  field  in  the  bore  of  the  field  winding 
is  shown  in  Figure  9  along  with  the  superconductor  short 
sample  test  curve.   Shield  structure  is  assumed  to  be 
stainless  steel  and  low  values  of  working  stresses  are 
chosen  due  to  the  fact  that  heat  dissipation  will  degrade 
material  properties.    The  120  degree  field  winding  angle 
has  an  advantage  of  eliminating  third  harmonic  component  of 
field  thus  preventing  flow  of  third  harmonic  currents  in 
delta  connected  machines.   In  addition  to  the  input  para- 
meters given  in  Table  3,  field  winding  inside  radius  of 
2.76  inches  and  outside  radius  of  5.0  inches  were  chosen, 
and  the  iron  shield  was  assumed  to  be  on  the  armature  end 
turns  (t  =  0) . 

Although  the  program  could  be  run  several  times  with 
■  different  field  winding  radii  and  input  parameters  to  get 
absolute  minimum  weight  and  volume  machine,  this  process 
was  omitted  here  and  only  a  design  for  a  given  set  of  input 
parameters  defined  in  Table  3,  was  obtained. 
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TABLE  3 

DESIGN  CONSTANTS 

DESIGN  HORSEPOWER=    20450.25 

MVA  RATING=    15-25 

POWER  FACTOR=1.00 

DESIGN  RPM=3600.0 

RATED  FIELD  CURRENT  DENSITY=1.15    10  8  A/M2 

RATED  ARMATURE  CURRENT  DENSITY=3.00    10  6  A/M2 

MAXIMUM  FLUX  DENSITY  IN  FIELD  WINDING=5.00     TESLAS 

FIELD  WINDING  ANGLE=120.00    DEG. 

ARMATURE  WINDING  ANGLE=60.00    DEG. 

ARMATURE  SPACE  FACTOR=0.27 

ARMATURE  FILAMENT  DIAMETER=0 . 04  0    IN. 

YOUNG  MODULUS  OF  SHIELD  SUPPORT  MATERIAL=29 • 0    10  6  PSI 

POISSON  RATIO  OF  SHIELD  SUPPORT  MATERIAL=0 . 30  0 

ELECTRICAL  CONDUCTIVITY  OF  SHIELD=6.0    10  7  MHO/M 

DENSITY  OF  SHIELD  SUPPORT  MTERIAL=0 . 290    LBS/IN3 

WORKING  STRESS  IN  TENSION=50000 . 0    PSI. 

WORKING  STRESS  IN  COMPRESSION=50000 . 0    PSI. 

WORKING  STRESS  IN  SHEAR=30000 . 0    PSI, 

CIRCUIT  BREAKER  OPENNING  TIME=10.0    CYCLES 

RATED  FIELD  CURRENT=680 . 0    AMP. 

ALLOWABLE  FIELD  CURRENT  RISE=420.0    AMP. 
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2B   Effect  of  Primary  Shield  Location 

Figure  10  shows  the  effect  of  primary  shield  location 
where  the  induced  currents  flov/  under  fault  conditions, 
for  different  air  gap  lengths.   Calculations  are  done  for 
the  same  set  of  input  parameters  as  defined  in  Table  3«   The 
case  where  the  primary  shield  is  inside  the  shield  struct- 
ure, produces  less  radial  stresses  than  the  case  where 
primary  shield  is  outside  the  shield  structure,  under 
three  phase  short  circuit  from  load,  thus  may  require  less 
shield  structure  and  airgap  length  and  therefore  less 
machine  volume.   Two  computer  outputs  given  in  Tables  4 
and  5,  show  this  effect.   As  can  be  seen,  the  results  are 
not  very  different  from  each  other.   Although  the  case 
where  primary  shield  is  inside  shield  structure  gives 
smaller  machine  volume  and  length,  the  direct  exposure  of 
the  shield  structure  to  the  trapped  magnetic  fields  may 
induce  higher  heat  dissipation  thus  increasing  the 
refrigeration  power  required. 
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•  TABLE  4 

FIELD  WINDING  OUTSIDE  RADIUS=    5-000    IN. 

FIELD  WINDING  INSIDE  RADIUS=    2.760    IN. 

ARMATURE  INSIDE  RADIUS=    6.451    IN. 

ARMATURE  OUTSIDE  RADIUS=    9-354    IN. 

SHIELD  INSIDE  RADIUS=    12.257    IN. 

SHIELD  OUTSIDE  RADIUS=    19-887    IN. 

TOTAL  LENGTH=58.848    IN. 

STRAIGHT  SECTION  LENGTH=    27-237    IN. 

MACHINE  VOLUME=    42.313    FT   3 

WEIGHT  OF  ACTIVE  PARTS=    l64l8.903    LBS. 

PRIMARY  SHIELD  IS  OUTSIDE  THE  SHIELD  STRUCTURE 

AIR  GAP  LENGTH=    1.4514    IN. 

PRIMARY  SHIELD  OUTSIDE  RADIUS=    6.364    IN. 

SHIELD  STRUCTURE  OUTSIDE  RADIUS=    6.364'   IN. 

PRIMARY  SHIELD  THICKNESS=    0.0    IN. 

SHIELD  STRUCTURE  THICKNESS=    1.3410    IN. 

SHIELD  DEFLECTION  TOWARDS  FIELD  WINDING  UNDER  THREE  PHASE 

FAULT=    0.0220    IN. 
SHIELD  DEFLECTION  TOWARDS  ARMATURE  WINDING  UNDER  THREE 

PHASE  FAULT=-0.0205    IN. 
SYNCHRONOUSREACTANCE  P. U. =0.2533 
TRANSIENT  REACTANCE  P. U. =0.1875 
SUBTRANSIENT  REACTANCE  P. U. =0.128 3 
MACHINE  LOSSES  PERCENT  OF  MACHINE  RATING=0.751 
MAX.  RADIAL  STRESS  UNDER  THREE  PHASE  FAULT=    3482.8    PSI 
MAXIMUM  TORQUE  UNDER  LINE  TO  LINE  SHORT  CIRCUIT  FROM  LOAD 

P.U.=    10.303 
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•  TABLE  5 

FIELD  WINDING  OUTSIDE  RADIUS=    5-000    IN. 

FIELD  WINDING  INSIDE  RADIUS=    Z.jGO         IN. 

ARMATURE  INSIDE  RADIUS=    6.306    IN. 

ARMATURE  OUTSIDE  RADIUS=    9.282    IN. 

SHIELD  INSIDE  RADIUS=    12.257    IN. 

SHIELD  OUTSIDE  RADIUS^    19- 887    IN. 

TOTAL  LENGTH=    57-950    IN. 

STRAIGHT  SECTION  LENGTH=    26.773    IN. 

MACHINE  VOLUME=    41.667    FT   3 

WEIGHT  OF  ACTIVE  PARTS=    I6169.683    LBS. 

PRIMARY  SHIELD  IS  INSIDE  THE  SHIELD  STRUCTURE 

AIR  GAP  LENGTH=    I.3O63    IN. 

PRIMARY  SHIELD  OUTSIDE  RADIUS=    5-026    IN. 

SHIELD  STRUCTURE  OUTSIDE  RADIUS=    6.29^1    IN. 

PRIMARY  SHIELD  THICKNESS=    0.0    IN. 

SHIELD  STRUCTURE  THICKNESS=   1.2680    IN. 

SHIELD  DEFLECTION  TOWARDS  FIELD  WINDING  UNDER  THREE  PHASE 

FAULT=    0.0245    IN. 
SHIELD  DEFLECTION  TOWARDS  ARMATURE  WINDING  UNDER  THREE 

PHASE  FAULT=-0.0229    IN. 
SYNCHRONOUS  REACTANCE  P. U. =0.2519 
TRANSIENT  REACTANCE  P.U.=0.l84  9 
SUBTRANSIENT  REACTANCE  P. U. =0.1657 
MACHINE  LOSSES  PERCENT  OF  MACHINE  RATING=0.749 
MAX.  RADIAL  STRESS  UNDER  THREE  PHASE  FAULT=    3377-9    PSI 
MAXIMUM  TORQUE  UNDER  LINE  TO  LINE  SHORT  CIRCUIT  FROM  LOAD 

P. U. =8. 061 
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2C   Comparison  Be-cween  20,000  HP  Conventional  Generator  and  a 
Superconducting  Generator 

The  results  obtained  In  previous  section  are  compared 

13 
with  a  20,000  HP  conventional  propulsion  generator   . 

Figure  11  shows  such  a  comparison  in  machine  dimensions. 
An  additional  15,000  lbs  are  added  to  the  weight  of  the 
superconducting  generator  in  order  to  take  account  of  the 
machine  casing,  bearings  excitation  system,  etc.   Table  6 
lists  some  of  the  information  available  about  the  conven- 
tional generator.   Unavailability  of  additional  information 
limits  comparison  in  machine  dimensions  and  weight.   As 
far  as  refrigeration  weight  and  volume  are  concerned,  there 
will  not  be  anything  to  say  unless  rotor  is  structurally 
and  mechanically  designed.   In  addition  refrigeration 
system  may  be  integrated  to  provide  for  other  superconduct- 
ing electric  machinery. 


TABLE  6 

20,000  HP  Conventional  Generator  Characteristics 

Maximum  Continuous        2  Minute 
Rating  Capacity 

KVA  15,250 

RPM  3,600  900 

Volts  4,700  823 

Power  Factor  1.0 

Frequency,  Hz  60 

Amperes  300^  Max 
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Chapter  3 
CONCLUSIONS  AND  RECOMMENDATIONS 

One  major  conclusion  that  can  be  drawn  from  the  previous 
results  is  that,  superconducting  machines  can  be  built  to 
withstand  fault  loadings  and  still  have  low  weight  and 
volume . 

Since  most  of  the  weight  comes  from  the  iron  shield, 
replacement  of  it  by  a  conductive  image  shield  may  lead 
to  additional  weight  savings,  although  image  shield  may 
introduce  additional  losses. 

Replacement  of  copper  by  aluminum  as  a  stabilizer  will 
also  reduce  field  winding  weight  if  a  suitable  fabrication 
technique  can  be  found. 

Conversion  of  the  computer  program  for  many  pole 
machines  will  increase  its  usefulness  and  provide  informa- 
tion about  low  speed  (100-200RPM)  synchronous  propulsion 
motors . 

Although  a  simple  technique  could  be  employed  which 
changes  input  parameters  systematically  and  gives  direction 
towards  minimum  weight  or  minimum  volume  machine,  it  is  the 
belief  of  the  author  that  such  technique  would  require  a 
great  amount  of  computer  time  and  money.   Furthermore, 
minimum  weight  or  volume  machine  may  not  be  best  as  far  as 
rotor  mechanical  design  and  refrigeration  input  are 
concerned. 

Perhaps,  the  most  questionable  assumption  is  that  the 
weight  densities  for  individual  machine  parts  are  uniform, 
especially  for  the  rotor  structure.   Correct  estimation  of 
rotor  weight  can  be  done  only  by  structural  and  mechanical 
design.   Critical  speeds  of  the  rotor  may  not  be  acceptable 
and  in  that  case  rotor  structure  or  length  should  be 
modified.   Rotor  lenght  should  be  increased  as  much  as 
possible  in  order  to  increase  thermal  resistance  and  thereby 
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reduce  heat  leak  into  cold  region.   Critical  speeds  will, 
however,  limit  the  rotor  length.   Inclusion  of  above  consider- 
ations will  require  longer  computation  time  but  may  give 
more  useful  information  about  superconducting  machines. 
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Appendix  A 
SHIELD  DEFLECTION  AND  STRESSES  UNDER  FAULT  CONDITIONS 

It  has  been  assumed  that  only  the  shield  structure  will 
be  effective  for  carrying  fault  loads.   Figure  A-1  shows 
radial  loads  imposed  upon  simply  supported  shield  structure 
under  three  phase  fault,  with  the  addition  of  centrifugal 
loading.   It  is  easier  to  separate  these  two  loadings  and 
superimpose  them  after  each  solution  is  obtained  separately. 

Reference  19  gives  three  differential  equations  for 
deflections : 


S^u   1-v  32 u  .  1+v   8^u    V  aw  _  p) 
9  x2    2a^  fe^   2a   8x39  ~  a  8x 


(A-la) 


1+v  ^u        1-v  8^u    1   9^u  _  1   gv 


2a   8x86    2   9x2    a^  90 


h' 


a''  99    12a' 


9^w 


+   33W 


9x^96    a^9 


+  h\ 
12a^ 


1-v  9^0  ^  __92u 


9x 


a^99 


=  0   (A-lb) 


V  ^+  ^^ 


9x   a9i 


w 
a 


hi 
12 


^  9x'*    a  9x^96"    a^9e'^ 


hi 
"12 


2-v   9'u 


1   9^u^ 


a   9x^90   a^  99 


^-^Eh =  0 


(A-lc) 


Simply  supported  boundary  conditions  imply  that. 


at  X  =  0  and  x  =  L 


u=0,w=0,  N^=0,  M  =0 
'        '   X     '   X 


These  boundary  conditions  can  be  satisfied  if  one  chooses 
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+  IJ2. 
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Figure  A-1 

Three  Phase  Fault  Radial  and  Centrifugal  Lc  -ling  on  Shield 

Structure 
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u  =  cos-^r   ^  A   cos  nQ  (A-2a) 

L    n 


u  =  sin^  E  B   sin  n0  (A-2b) 

L    n  ^     ' 


w  =  sin^  Z  C   eos  ne  (A-2c) 

L    n 

It  was  shown  in  Chapter  2-F  that  the  radial  stresses 
depend  on  the  angular  position  along  the  circumference. 
They  are  given  by  equation  21. 

2    ^o 
Po  =  a  =   a     cos  9  =  ^5—  (1  +  cos  26)  x  sin— :r   ,„  .,  v 
3    r    o         2  L   (A-3) 


When   equations    (A-2)    and    (A-3)    are    combined   with    (A-1), 
and  the    last    two   terms   neglected   in   equation    (A-lb)    because 
they   will   be    small    for   thin   shells    compared  to   the   other 
terms, 

A    [2Trn(l-v)X^n^]-B    [  (1+v)  X7Tn]+C    [2vA7t]  =  0  (A-4a) 

A    [3(l+v)7TnX]-B^[3(l-v)TT2+6A2n2]+C^[6A2n]  =  0  (A-4b) 

A^[v7t3X]-B^   nX2r3+Ti2[(2-v)7T2  +  x2n2]    +C^[3X2+n' (Tr^  +  X^n^  )  ^] 

=  ^P„  (A-4c) 


where  X  =  — 


^  =  2L 


h_ 

3 


D  =  -^^ 


12(l-v2) 
a. 


and  for  n=0  and  2,  P  ~^2~?~'    ^^^  result  which  comes  directly 
from  equation  (A-3) 
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The  coefficients  obtained  by  solution  of  equation  (A-4) 
are  too  complicated  to  write  explicitly.   The  results  are 
summarized  in  Table  A-1.   The  stresses  can  be  obtained  by. 


^x   izf^  '^^x'^^^e'^^^x'^^'^a^^  (A-5a) 


^6  =  irfa-  [sQ+ve^-z(XQ+vX^)]  (A-5b) 


^xe  =  mwy  ^'^-^^^xe)  ^^-5c) 


where 

9u 


^x    8x 


En  ^  1  8_u    w 
a  80  ~  a 


1  ,1  ?u  .  9^ 


w 


(A-6a) 


(A-6b) 


.    8u  ,  1  9u  f^    c    \ 


\   =  ^  (A-6d) 


V    1  3   ,u  ,  1  3w  ■  ,.    £    s 

^=9  =  a  36-  (a  +  a  30-'  '*-^^' 


\e  =  I  ^IW^   33ffF)  ('^-Sf ) 


and  z  is  the  distance  from  the  mid-surface  (positive 
towards  the  axis  of  the  shell).   Equation  (A-6)  can  be 
combined  with  equation  (A-5)  to  give  stresses  in  axial  and 
circumferential  directions. 
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TABLE   A-1 
Solution   of   equation  A-4 


A      = 
o 


^  C 

IT  O 


Bo    =    0 

^o   "   2D[3A^(l-v'^)+n'7T'^] 


A^    =    -B^ 


^12X^l±y)__    3(1.^)^2.24,.^ 


[6(1+v)7tX-^[7t2+(i-v)2aF 


■p       —    ; 

2  2D 


6(1+v)ttA-^^[7t2+(1-v)2X2] 

V7T 


/ 


1  +  V 


(TT2  +  i|X2)]-2A2[3  +  n2[(2-V>TT2+i|X2]] 


6(l+v)7rA- 


[3A^+n' 

12A 


VTT 


[tt^+(1-v)2A^] 


"^^^'(^^^)      -3(1_v)7t2-24A^' 


BvttA- 


[27t2+(1-v)4a2]    .     [3X2+nMTr2+4A2)2]/[2vA7T] 


^2    = 


B2    .    [2(1  +  v)A7t]-A2    .    [27r2+(l-v)4A2] 

2vA7T 
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Since  the  shield  will  rotate  with  field  winding, 
centr:.fugal  stresses  will  be  induced  with  associated 
deflections.   The  centrifugal  loading  may  be  considered  as 
loading  due  to  uniform  internal  pressure  given  by: 

pw^.ah 

i~  (A-7) 


P  = 


Reference  19  gives  the  solution  for  a  circular  cylind- 
rical shell  with  a  uniform  internal  pressure  with  ends 
simply  supported: 


^  -  -  Md^ 


1-  7:: 7 T-TT-   .  sin3(x-Tr)  .  sinh3(x-7r) 

cos2o=  +  cosh2a        ^   2         ^^      2 


2cos°^ 


cosh   °^ 


cos2<3:   +   cosh2°: 


cos6(x-|) 


cosh3  (x-p-) 


(A-8) 


where 


3  =     ^  ^^"h^^ 


and 


«  =  3L 
2 


Equation  (A-8)  is  the  deflection  in  vertical  direction 
Prom  symmetry,  the  component  v of  the  displacement  in  the 
circumferential  direction  vanishes.   The  expressions  for 
the  strain  components  then  become 


du 
dx 


(A-9a) 


— 
^e   "a 


(A-9b) 


and  it  can  be  shown  that 

du     w 
d7  =  ^  a 


19 


(A-10) 
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Equations  (A-5),  CA-S ) , (A-7 ) ,  (A-8),  (A-9),  (A-10) 
and  Table  A-1  can  be  combined  to  get  total  deflections 
and  stresses  along  the  shell  at  any  point  on  circumference. 
The  component  w  of  the  displacement  will  be  maximum  at 
midlength  (x=L/2)  when  9  =  0  and  9=7t/2  denoted  by  w,  and  Wp 
respectively.   These  two  deflections  will  determine  the 
required  airgap  along  with  shield  thickness.   If  the 
material  properties  are  chosen,  it  is  possible  to  calculate 
stresses  at  these  points,  top  and  bottom  surfaces.   Since 
the  shear  stresses  will  be  zero,  these  are  actually  prin- 
cipal stresses.   Maximum  shear  stresses  can  be  obtained  by. 


a  -a, 
X   ( 


xe 


(A-11) 


Under  line-to-line  short  circuit  from  load,  large  fault 

torques  are  applied  on  the  shield.   The  magnitude  of  these 

2^ 
torques  is  given  by   .   The  shield  structure  must  be  strong 

enough  to  carry  these  stresses.   For  the  shear  stress  in  the 

23 
she-11,  according  to  the  maximum  shear  theory 


(•'xe^   =  2 

max 


^^-xe^' 


1/2 


(A-12) 


in  this  case   6  will  be  assumed  to  consist  of  centrifugal 
stresses  only,  that  is. 


2   2 

pw  r 


(A-13) 


and 


xe 


7t(R 


2R   T 
so 


so  -  si  ) 


(A-14) 


Equations  (A-13),  (A-l4),  and  (A-12)  may  be  combined  to  give 
maximum  shear  stresses. 
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Application  of  this  fault  torque  may  cause  buckling  of 
shield.   Reference  20  gives  criteria  for  critical  buckling 
stresses . 


2/^o^.^   ,.2n-5/8^,^,„,5/^,_  .r  ^l/2 


•   T^^  =  0.85  (7tV12)(1-v^)  ^'^E(h/a)-''"(a/L) 
for  10(h/a)-^^^<L/a<3(a/h)^'^^ 


for-L/a>3(a/h) 
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(A-15a) 


T^^  =  0.272  (l-v^)-3/^E(h/a)3/2  (A-15b) 

1/2 


Appendix  B 

COMPUTER  PROGRAM  LISTING  AND  DESCRIPTION 

The  first  harmonic  field  in  the  field  winding  for  a 
two  pole  (p=l)  machine  is  given  by  (Table  1). 


B 


2y^  J^  sin 

0    I 


rfw 


wfe 


sin(e-(^) 


3tt 


-2- 

r 

3 
+3 

r 

hfoj 

-1 


r 

-1 

(i-y^) 


(A-16) 


The  maximum  flux  density  in  the  r  direction  can  be 

found  by  taking  the  derivative  of  B  „  with  respect  to  r 
^  rfw         ^ 

and  setting  it  equal  to  zero,  which  gives  a  simple  result 
for  a  two  pole  machine 


r  =  R 


fi 


(A-17) 


and  the  maximum  flux  density  in  the  field  winding  is 

e. 


B 


2vi   J^  sin 
^0   f 


wfe 


R 


fo 


max 


Stt 


-3y+3+ 


^^fo^2 


R 


(1-y^) 


(A-18) 


This  linear  relation  between  B    and  J^  defines  the 

max      f 

magnet  line  which  was  shown  in  Pigs.  8  and  ?•   If  an  oper- 
ating point  is  chosen  (given  B  ^   and  J^)  and  field  winding 

max      I 

inside  and  outside  radiuses  are  specified,  equation  (A-l8) 

can  be  solved  for  R  ,  shield  inside  radius.   Although 

fields  generated  by  the  armature  phases  try  to  oppose  the 

fields  generated  by  field  winding  thus  reducing  B    ,  there 

max 

will  be  additional  field  concentration  due  to  armature  end 

turns  and  higher  space  harmonics.   Thus,  calculating  B 

by  equation  (A-l8)  is  a  conservative  assumption. 
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After  the  location  of  the  iron  shield  is  fixed,  the 
first  thing  to  be  checked  is  the  saturation  limitations. 
Fig.  A-2  shows  the  flux  linkages  through  the  iron  shield. 
The  total  flux  carried  by  the  shield,  for  a  two  pole 
machine 


•tt/ 


$,  ,  ^r^"    B  .   (RJ  I      R   d 


rfo    s^  *s   s 


(A-19) 


where  the  first  harmonic  component  of  B  ^   comes  from 

rfo 

Table  1  and  is  given  by 

'0.   ^ 
2J^  sin 
B 


wfe 


cos6 


rfo 


37T 


R 


fo 


(1-y^) 


1+ 


r 


^. 


(A-20) 


When  equation  (A-20)  is  evaluated  at  r  =  R  and  substituted 
into  (A-19),  integration  gives 


Va  =  ^%  "^f  ^^"^ 


wfe 


3tt 


^^5,(1-^3) 


(A-21) 


The  limitations  on  saturation  require  that 


^A-A  -<  ^sat  ^^os-f^s^  ^s 


(A-22) 


When  equations  (A-21)  and  (A-22)  are  combined,  the  minimum 
shield  outside  radius,  R^^  is 

'    OS 


R 


4y^  J^  Sin 


wfe 


OS 


Btt  B 


sat 


R^  (1-y^)  ^  % 

s 


(A-23) 


where  B   ,  is  the  saturation  flux  density  of  the  iron 
sac  P 

shield.   B   ,  =  1.5  Weber/m  was  used  in  the  program, 
s  a  T^ 

Local  saturation  of  the  iron  shield  requires  that 


Brfo(Rs)  = 


4y   J^  sin 
0   f 

3^ 


fe„fel 

2 

D 

[«fol 

"s 

["^  ] 

(1-y^)  ^  B 


sat    (A-24) 
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Figure  A-2 
Flux  Pattern  in  Iron  Shield 


Algebraic  manipulation  of  equation  (A-23)  and  (A-24) 


gives 


R   -  R  ^  R 
OS     s     s 


(A-25) 


This  is  to  say  that  shield  thickness  must  be  less  than 

or  equal  to  shield  inside  radius.   Since  R   is  fixed  from 

s 

equation  (A-I8),  calculated  R   may  not  satisfy  this 

o  s 

condition.   In  this  case,  program  rejects  this  particular 
design  and  prints  a  message  saying  "NO  SOLUTION".   Same 
message  is  also  printed  when  equation  (A-I8)  gives  R  less 


than  R 


fo* 


The  maximum  airgap  length  can  be  found  by  assuming 
that  armature  inside  and  outside  radiuses  coincide,  that  is 


A    =  R  -t-R^ 
max    s    fo 


(A-26) 


where  t  is  the  radial  distance  between  shield  inside  radius 

and  armature  end  turns.   In  the  program,  this  A    is 

^   *^    '       max 

64 


divided  into  n  steps  (n=50  in  this  program)  and  in  each 

iteration  airgap  length  A  is  incremented  by  A   /50. 

max 

Given   A,  armature  inside  radius  R  . 
'  ai 

\i   =  Rf„  +  A  CA-27) 

and   ai''mature   outside   radius,    R 

'   ao 

Ro^  =  (R^  +  R  .  -  t)/2.0         (A-28) 
ao     s    ai 

The  program  then  proceeds  to  sum  up  the  infinite 
series  given  in  Table  A-2  of  Chapter  2  for  L   -  L  ,  .   Then 
synchronous  reactance  is  normalized  to  internal  voltage 
and  X   as  given  by  equation(6 )is  calculated.   Table  A-2 
gives  the  effective  lengths  used  in  this  study   .   Equation 
7  is  then  utilized  to  give  an  expression  for  the  straight 
section  length,  which  results  in  a  second  degree  equation 
for  the  straight  section  length.   This  equation  is  too 
complicated  to  write  explicitly  and  will  be  omitted  here. 
Total  machine  length,  is  given  by 


I      =    i    +   2A£  (A-29) 

After  all  machine  dimensions  are  calculated,  the 
program  then  proceeds  to  calculate  machine  reactances. 
Synchronous  reactance  x.  given  by( 9) , transient  reactance 
X . '  given  by(8)can  be  calculated  by  summing  the  series.   It 
should  be  pointed  out  that  up  to  15  harmonics  were  used. 
The  user  is  free  to  choose  number  of  harmonics  to  be  summed 

To  calculate  subtransient  reactance  x,"  ,  the  primary 
shield  radius  has  to  be  specified.   Since  for  a  given  air- 
gap  length,  the  range  of  R.  is  known  and  this  can  be  done 
easily  by  incrementing  R,  in  each  iteration.   For  a  prev- 
iously calculated  airgap  length  A,  R  is 
(m=50  used  in  this  study)    R^  =  R^^  ^  A  (A-30) 
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TABLE  A-? 

Effective  Lengths 

Armature  end  turn  length  AJl   =  R  .  +  R 

^  ai    ao 

Effective  lengths:   armature        i  o  .  ao 

-,  ^    .     ■,      ,  a=X'  +  Ax, 

self  inductance 

Mutual  inductance,  field-to- 
armature 

Mutual  inductance,  field-to- 
damper 

Eddy-current  loss 

Conduction  loss 

Field  winding  end  turn  length 


z     = 
m 

i 

^fd   = 

Z 

^d   = 

Z 

^c    = 

Z+    2A£ 

A£^   = 

^fo    ^   ^fi 
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Then  x . "  can  be  calculated  from(lO). 

Next  step  Is  to  calculate  required  primary  shield 
thickness,  to  give  time  constant  such  that,  field  current 
rise  will  not  be  higher  than  the  specified  value  (SUB- 
ROUTINE TIC).   Since  the  maximum  field  current  rise  will 
be  from  2  5 


=  1 


max 


fo 


-X 


1  + 


(A-31) 


which  depends  on  x ^-x. '  .   If  the  maximum  rise  is  less  than 

X  ' 
^d 

the  specified  rise,  then  the  required  shield  thickness  or 
time  constant  is  set  to  zero.   This  is  to  say  that  the 
primary  shield  will  be  a  very  thin  highly  conductive  layer 
where  induced  currents  flow.   In  this  study,  it  was  assumed 
that  the  time  constant  of  shield  structure  is  zero. 

Having  specified  primary  shield  thickness,  the  radial 
stresses  a   under  three-phase  fault  from  load  is  calculated 

0 

(SUBROUTINE  STRES)  by(21). 

The  most  important  section  of  the  program  is  SUBROUTINE 
SHIEL  which  designs  the  shield  structure  according  to  the 
given  working  stresses.   The  designed  shield  must  satisfy 
the  following  conditions: 


1) 


W-,  +  h,  +  h  +  w„  =  A 
1    h        2 


where  h,  is  the  primary  shield  thickness,  h  is  the  shield 
structure  thickness  and  w-,  and  Wp  are  the  deflections  in 
two  directions. 


2) 


\    -   \-   ^/2    =  ^ 


\  +  h/2  =  a 
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-if  primary  shield  is 
outside  the 
shield  struct- 
ure 

-if  primary 
shield  is  inside 
shield  structure 


3)  S^  -^  "sic 
St  «  "sit 
Ss  *  "sis 


where  S  ,  S,  ,  S   are  induced  stresses  in  compression,  tension 
c    t '   s 

and  shear  respectively  and  W  .  ,  W  .,  ,  W  .   are  the  working 

^         "^  SIC    sit     SIS  '=' 

stresses  in  compression,  tension,  and  shear  respectively. 


x6        sis     stresses  induced  under  line- 
max 

to-line  short  circuit  from 
load 

"^cr  "^  ^-^xe^ 

max 


A  total  of  8  principal  stresses  are  calculated  after 
W-,  and  W^  are  found  by  equation  (A-4).   These  are  at  0=0 
and  6='n'/2,  on  the  top  and  bottom  surface.   All  of  the  stres- 
ses mentioned  above  are  at  the  midlength  where  maximum 
stresses  occur.   Shield  length  was  assumed  to  cover  straight 
section  length  plus  field  winding  end  turns.   After  checking 
each  of  them  by  the  requirements  specified  above,  the 
progrm  proceeds  to  check  on  requirement  4.   If,  in  the 
course  of  these  calculations,  any  one  of  the  requirements 
is  not  satisfied,  computations  are  repeated  by  incrementing 
shield  structure  thickness,  until  specified  airgap  is  filled 
with  shield  structure.   Figure  A-3  shows  the  region  to  be 
scanned  by  the  program.   Since  the  smallest  airgap  length 
will  give  minimum  weight  or  volume  machine,  the  solution 
will  return  to  the  main  program  if  all  four  of  the  require- 
ments are  satisfied.   There  may  be  other  solutions  at 
bigger  airgap  lengths  which  we  will  not  investigate. 
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Machine  volume  is  calculated  by 


V  =  7T  R  ^ .     i. 
OS      t 


and  machine  weight  by 


where 


-  armature  density        =  ^500  Kg/m-^ 


Pa 

3 

p   -  iron  shield  density     =  8000  Kg/m 
s 

3 

p„  -  rotor  density  =  ^500  Kg/m 

3 
p   -  end  turns  density       =  4000  Kg/m 

p   -  shield  structure  density=depends  on  program 
^^  input 

These  average  densities  are  used  in  the  program. 

For  selected  values  of  n,  m  and  increment  of  h,  the  shield 

structure  thickness  solution  may  not  converge  or  may 

converge  at  a  longer  airgap  length,  if  these  values  are 

small.   Smaller  step  sizes  will  increase  computation  time 

and  give  true  convergence,  so  some  sort  of  a  comprimise 

should  be  done.   If  the  solution  does  not  converge  then 

program  prints  a  message  saying  "NO  CONVERGENCE". 
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ct 


^k'%- 


An,ax/n 

i 


_  X  Amax/nm  \ 


I  \ 


'    A  =  4A,ax/n 


A  =  3A«ax/n 


\   A  =  2A„ax/n 


\  A=  A^.Vn 


max' 


■v    A=0 


Figure  A-3 
Region  Scanned  by  the  Program 

Main  program  and  subroutines  are  listed  in  Tables  A-3, 
A-4j  A-5,  and  A-6.   Table  A-6  gives  the  sample  program 
input  which  is  described  below. 


Card  1   Format  3F10.2 


Horsepower  Rating 
Power  Factor 
Design  RPM 


col.  1-10 
col. 11-20 
col. 21-30 


Card  2   Format  2F10.2 


Field  winding  outside  radius,     col.  1-10 
inches 


Field  winding  inside  radius, 

inches 
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col. 11-20 


Card  3  .  Format  3F10.2 

Maximum  flux  density  in  field     col.   1-10 
winding,  Tesla 

Ratednfield  current  density       col.  11-20 
X  10   ,  amp/m 

Field  winding  angle,  0,,^  col.  21-30 

T  WI 

degrees 

Card  4   Format  4P10.2 

Rated -armature  current  density    col.   1-10 
X  10~°,  amp/m2 

Armature  winding  angle,  9   ,      col.  11-20 

-,  wa 

degrees 

Armature  space  factor  col.  21-30 

Armature  filament  diameter,       col.  31-^0 
inches 

Card  5   Format  7F10.2 

Shield  structure  material         col.   1-10 
young  modulus  x  10~",  psi 

Shield  structure  material         col.  11-20 
poisson  ratio 

Primary  shield  conductivity       col.  21-30 
X  10""^,  mho/m 

Shield  structure  density,  lbs/in^  col.  31-^0 

Working  stresses  in  tension,  psi   col.  41-50 

Working  stresses  in  compression,   col.  51-60 
psi 

Working  stresses  in  shear,  psi    col.  61-70 

Card  6   Format  3F10.2 

Circuit  breaker  opening  time,      col.   1-10 
cycles 

Rated  field  current,  amp  col.  11-20 

Allowable  field  current  rise, amp   col.  21-30 

71 


"Card  7 

Format  II,  if  1  is  punched  in  column  one,  program 
makes  calculation  for  the  case  where  primary 
'shield  is  under  shield  structure;  and,  if  2  is 
punched,  calculations  are  performed  for  the  case 
where  primary  shield  is  on  the  shield  structure. 

Format  F10.2,  distance  between  shield  inside 
radius  and  armature  end  turns,  inches    col.  11-20 
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